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Abstract 

Oh" 

Inhomogeneous fluid flows which become supersonic are known to produce acoustic 
analogs of ergoregions and horizons. This leads to Hawking-like radiation of phonons 
with a temperature essentially given by the gradient of the velocity at the horizon. We 
find such acoustic dumb holes in charged conformal fluids and use the fluid-gravity cor- 
■ respondence to construct dual gravity solutions. A class of quasinormal modes around 

CN ■ these gravitational backgrounds perceive a horizon. Upon quantization, this implies a 

thermal spectrum for these modes. 
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1 Introduction 



In recent years, gauge-string duality [U E] has been useful in exploring properties of strongly 
coupled field theories in regimes where their duals may be truncated to classical gravity. In 
particular, application of gauge-string duality to the hydrodynamic regime of these field the- 
ories has led to a "fluid-gravity correspondence" [3]- [7]. Properties of solutions of classical 
gravity then lead to predictions for interesting properties of the dual fluid, the most celebrated 
example being the ratio of shear viscosity to the entropy density of a conformal fluid [3J. 

In this note we use the fluid-gravity correspondence in the opposite direction. We use 
properties of supersonic fluid flows to predict interesting properties of fluctuations around a 
class of deformed black brane spacetimes in asymptotically AdS spacetimes. These spacetimes 
are duals of inhomogeneous flows of conformal fluids where the fluid velocity exceeds the speed 
of sound in some region. Unruh [8J showed that such flows lead to the formation of an "acoustic 
ergoregion" and, under suitable conditions, to an "acoustic horizon" . The same physics which 
leads to Hawking radiation from black holes in General Relativity now leads to a Hawking-like 
radiation of quantized sound waves (or phonons) with a thermal spectrum, the temperature 
being proportional to the gradient of the velocity field at the acoustic horizon [TU| [TTJ . Even 
when an acoustic horizon is not present, the presence of an ergoregion leads to characteristic 
properties like superradiance [13] . Fluid configurations with such acoustic horizons have been 
termed "dumb holes", and have been proposed as possible experimentally realizable systems 
for testing the physics of Hawking radiation in the laboratory [13] . 

We will show that the gravity duals of such supersonic flows are non-static black holes. 
The duals of sound waves are then certain quasinormal modes around such black holes, and 
it follows from the fluid-gravity correspondence that at the quantum level one should find 
a Hawking-like radiation of these modes with an approximately thermal spectrum [3]. This 
Hawking-like radiation is distinct from the usual Hawking radiation associated with the black 
hole horizon, and would be present even when the background black hole is extremal and 
hence at zero temperature. The temperature of this quasinormal mode radiation depends on 
the properties of a "quasinormal mode horizon", which is an extension into the bulk of the 
acoustic horizon of the boundary fluid. 

It should be emphasized that this phenomenon could have been found purely in General 
Relativity (or its supergravity extensions relevant to our considerations) by studying the fluc- 
tuation problem around these non-static black holes. However, without the fluid gravity corre- 
spondence and knowledge of acoustic Hawking radiation, there would not have been an obvious 
motivation to look for quasinormal mode horizons in non-static black brane backgrounds. 

While we believe that the phenomenon of Hawking radiation or super-radiance of quasi- 
normal modes is quite general, it turns out to be rather difficult to come up with examples 
within a controlled approximation scheme. The simplest and perhaps most interesting back- 
ground where such a phenomenon could be present is a Kerr black hole in asymptotically AdS§ 
spacetime. The dual of such a background is a rotating conformal fluid on S 3 [T7J [T_B] . There 
is a regime of parameters of the black hole geometry for which the dual rotating fluid has 
supersonic velocities in a band around the equator of the S 3 , thus producing an ergoregion for 
sound modes. The physics of sound waves around such a rotating fluid background would have 
a dual description in terms of quasinormal modes of gravitational perturbations around the 
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Kerr black hole in AdS§. However, this flow has vorticity, and the existence of acoustic Hawk- 
ing radiation has been demonstrated mostly for irrotational flow. In the presence of nonzero 
vorticity, the sound modes get mixed up badly with other modes and analysis becomes difficult 

The situation simplifies, however, if the flow is irrotational. As shown in |8J for non- 
relativistic perfect fluids, and extended to relativistic perfect fluids in [20], the velocity poten- 
tial then obeys the wave equation for a minimally coupled massless scalar field propagating on 
a curved background, the metric of which is determined by the underlying flow. The mathe- 
matical problem of quantizing sound waves or phonons around such a flow is then quite similar 
to that of quantizing a massless scalar field in an ordinary black hole background. This implies 
the existence of an acoustic analog of Hawking radiation. 

Known examples of supersonic flows of perfect fluids often lead to infinite "acoustic surface 
gravity" (which is proportional to the gradient of the velocity at the acoustic horizon). The 
presence of viscosity usually regulates this divergence and renders it finite [IB]. The incorpo- 
ration of viscosity, however, makes the analysis complicated. 

In this paper we find simple examples of acoustic horizons in ideal relativistic conformal 
fluids with finite acoustic Hawking temperature. The simplest example involves a fluid moving 
in a background spacetime of the form 

ds 2 B = -dt 2 + dz 2 + R{z) 2 {d6 2 + sin 2 6d<f) 2 ) (1) 

where R{z) is a slowly varying function which has the behavior R(z) — > \z\ as \z\ — > ooJ§ 
The fluid flow is steady and the only nonzero component of the fluid velocity is v z (z), with all 
derivatives bounded. Starting with v z = at z = — oo we will show that v z reaches the speed of 
sound - producing an acoustic horizon - at minima of R(z). If the function R(z) has only one 
minimum, e.g. R(z) = \J (z 2 + z 2 ), the assumption of a smooth solution of the fluid equations 
of motion implies that the fluid velocity continues to increase beyond the acoustic horizon 
until it reaches the speed of light at z = oo. However, if R(z) has multiple extrema, e.g. two 
minima separated by a maximum, we will find smooth flows where the fluid velocity reaches a 
maximum supersonic value at the maximum of R{z) and then decreases, then turns subsonic 
at the second minimum, and finally reaches zero at z = oo. Sound waves cannot escape to the 
asymptotic region z = — oo from beyond the first acoustic horizon (which is therefore like a 
black hole horizon), and cannot cross the second acoustic horizon (which behaves more like a 
white hole horizon) from the z = oo asymptotic region . 
We also study flows in a warped R 1 ' 1 x T 2 geometry 

ds 2 = -dt 2 + R{z) 2 {d6l + de 2 2 ) + dz 2 (2) 

and find very similar phenomena. 

The acoustic Hawking radiation that arises when the sound modes of these flows are quan- 
tized will be easiest to detect if the acoustic Hawking temperature T H is larger than the 
ambient temperature of the fluid, Tjj > T. For an uncharged conformal fluid, the only scale is 

5 The precise meaning of a slowly- varying R(z) is given in the discussion above E cm. (1871) . R{z) is some 
function in the metric (TT]), for which all invariants constructed out of the curvature and its derivatives are small 
compared to the basic scale in the problem. 
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the temperature, so that hydrodynamics is valid only when all derivatives are small compared 
to the temperature. However Tjj itself is proportional to the gradient of the velocity field, 
Th ~ ^-\ z =z- Thus it is not possible to have Hawking radiation at a temperature higher 
than the ambient temperature in an uncharged conformal fluid. We will therefore look for 
charged fluid solutions satisfying the rather stringent criterion of Th > T, in order to ensure 
detectability of Hawking radiation]^ Note that the solutions we will find all have well-defined 
uncharged limits and satisfy all criteria other than Th > T under the replacement of T — > T 
and q — > 0. 

We will consider conformal fluids with global charge density q and conjugate chemical 
potential fi. For such a fluid in equilibrium, the ambient temperature can be made to vanish 
provided the charge density is chosen properly. In analogy with the gravitational duals of such 
fluids considered below, we will call such a fluid "extremal" . We will consider isentropic flows of 
such a fluid where the charge density q, energy density e and the velocity vary slowly (in a sense 
defined precisely below) but whose variations can themselves be 0(1). In an isentropic flow, the 
entropy density per unit charge density is, however, a constant. For such flows q oc e 3 / 4 where 
e is the energy density. Then T oc is the only independent energy scale in the theory. T is 
in general a function of the chemical potential // and the temperature T. As shown below, the 
isentropic condition allows us to keep the local temperature T to be always much smaller than 
T, even though the other hydrodynamic quantities can change by 0(1). In the limit of very 
small temperature T oc //. One would expect that the hydrodynamic approximation is valid so 
long as all gradients are small compared to T. We will show that it is consistently possible to 
construct fluid flows described above with all gradients ^ and all curvature invariants much 
smaller than T, thus ensuring T T H ^> T. 

The spacetimes (CD) and (j2J) may each be regarded as the boundary of an asymptotically 
AdSs near-extremal charged black brane geometry which is deformed due to a nonzero bound- 
ary curvature. For generic spacetimes with arbitrary R(z), the boundary metric might not 
admit a smooth bulk dual. However it has been shown in [7] that for slow-enough variations 
in the boundary metric, the dual geometry is regular up to the bulk horizon and free of other 
singularities. Our solution will admit a smooth bulk dual as long as all invariants constructed 
out of the boundary curvature and its derivatives are much smaller than the radius of the 
outer horizon R + in AdS units. One should note that in spite of being slowly- varying, the de- 
formations can still be large. Very close to extremality, R + ~ T ~ fi, so that this condition 
is in fact the condition for validity of hydrodynamics in the boundary theory. In the second 
example, the boundary metric has two compact directions. This means that the nature of the 
dual geometry depends on the size of the compact directions compared to R + [23J. We will 
choose R(z) ^> R+ so that the dual is a near-extremal black brane rather than an AdS soliton 
with a small temperature. 

A fluid flow profile in the boundary theory is then described by a normalizable deformation 
of this bulk metric. We construct the deformed bulk metric using a derivative expansion, 
following [5], [21] and [22]. The straightforward derivative expansion breaks down in "tubes" 

6 We would like to thank the referee and also Dileep Jatkar for bringing our attention to the work of 
Weinfurtner et.al. [26] where analogue Hawking radiation at a very low temperature has been observed in a 
background that has a temperature several orders of magnitudes higher than the radiation. 
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of constant retarded time where the geometry becomes exactly extremal; therefore, we consider 
fluid flows where the local temperature is small but nonzero. We then consider the class of 
linearized fluctuations around this background geometry - quasinormal modes - which are dual 
to sound waves in the presence of the corresponding fluid flow. Note that while the deformations 
of the bulk metric due to a nontrivial R(z) and a nontrivial velocity profile v z (z) are typically 
large, the quasinormal mode amplitudes are small. 

The behavior of the quasinormal modes clearly shows that at leading order in the derivative 
expansion, the acoustic horizon of the fluid extends into the bulk in the following sense. Let 
r denote the radial coordinate in the AdS space and let z = z be the location of the acoustic 
horizon in the boundary flow. We find that for any value of r, these quasinormal modes suffer 
an infinite blue-shift as we approach z = z: modes which travel along the direction of the fluid 
flow are smooth at z — z, while the modes which travel in the direction opposite to the flow 
have rapid oscillations. Thus, in an eikonal approximation these quasinormal modes cannot 
cross the quasinormal mode horizon at z — z, which extends radially from the acoustic horizon 
into the bulk. 

Standard arguments imply that upon quantization^, one would find a thermal distribution 
of these quasinormal modes with a temperature Th, which is the gravity dual of the acoustic 
Hawking radiation in the fluid. Only these specific quasinormal modes perceive the quasi- 
normal mode horizon; other modes can cross it with ease. By the same token, the thermal 
distribution will be made up only of these quasinormal modes; it exists independently of (and 
at a different temperature from) the usual Hawking radiation associated with the event horizon 
of the background black brane. 

Our discussion is restricted to the lowest non-trivial order in the derivative expansion, 
which is consistent with the perfect fluid approximation. However we expect that the physical 
consequences should survive higher derivative corrections. Furthermore our discussion of fluc- 
tuations, both in the boundary fluid and in bulk gravity, is restricted to the linearized limit. 
We do not address the effect of nonlinear interactions of the sound waves and other modes. 

Admittedly, our setup is a bit contrived and is meant to provide a simple toy model in 
which this novel gravitational phenomenon can be studied in a controlled fashion. We expect, 
however, that the phenomenon is quite general and would be present in more interesting 
situations (e.g. the Kerr black hole mentioned above). 

The paper is organized as follows. In Section [21 we give a self-contained discussion of 
acoustic metrics and dumb holes for conformal relativistic fluids. In Section [3], we describe the 
bulk dual. In Section H] we discuss the regime of validity of our solutions. 



2 Acoustic metric for relativistic conformal fluid 

In this section we derive the equation governing the propagation of sound around gradient 
flows of a perfect relativistic conformal fluid. For such fluids, the pressure p and the energy 
density e are related by 

>=!■ (3) 

7 Quantization of bulk modes corresponds to l/N corrections in the SU(N) gauge theory on the boundary. 
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For a charged conformal fluid with charge density q, there is an additional equation of state 
e = e(s,q), or equivalently s = s(e,q), that relates the energy, entropy and charge densities. 
We will eventually be considering charged fluids with gravity duals, and will write down an 
explicit equation of state for such fluids in Section (J3J). 
The first law of thermodynamics reads 

de = T ds + (jl dq , (4) 

where T and \i are the intensive quantities temperature and chemical potential respectively 
and can be obtained from the equation of state by taking derivatives: 



j, de(s,q) 



ds 



_ de(s,q) 
[i - 



dq 



(5) 



For a homogeneous system, it follows from extensivity that all thermodynamic variables are 
related by a Gibbs-Duhem relationship, which using ([3]) may be written as 



e = Ts + jj,q . (6) 



In the following, we will define a quantity T with dimensions of energy by 

P= € - = cT\ (7) 

where c is a dimensionless constants depending on the underlying system. T will be a function 
of T and \x (or q), which reduces to T in the uncharged limit. Our fluids will also admit 
a zero-temperature, finite-// limit, close to which T is proportional to /a. T sets the energy 
scale of our conformal fluid, and will play an important role in defining limits in which our 
approximations are valid. 

The equations of motion of fluid dynamics are conservation of the energy momentum tensor 
and conservation of the currents associated with any conserved charges, including the conserved 
particle number: 

= 

and V M jf = 0. (8) 

The stress tensor T^ v depends on the 3 independent components of velocity v(x M ), the energy 
density, the pressure and their derivatives. The currents jjf additionally depend on the densities 
qi of the conserved charges. To leading order in the derivative expansion, 

T i» = p g ^ + ( e + p)u^u u = cT 4 (g^ + AuV) (9) 
J'f = (10) 

Here u M = (7,7V) and 7 = a . Conformal invariance implies that the stress tensor is 
traceless. In a general curved background there is a trace anomaly; however this is a higher 



3 For an uncharged conformal fluid, such a relation is trivial and s ~ es 
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order effect in the derivative expansion and we have ignored it. We have also ignored the 
viscous and diffusive terms - which are again higher order in the derivative expansion - and 
we therefore work in the perfect fluid limit. In addition we will also restrict attention to the 
case of a single charge of density q(x^). 

The parallel component of the equations of motion (jSJ), u u V = 0, leads to the conser- 
vation law 

v,(rV) = o. (ii) 

In the uncharged case T 3 is proportional to the entropy density of the fluid and the above 
equation is the conservation of the entropy current. The perpendicular component P^'V fl T fJ ' l/ = 
(where the projector P x = 5^ + u x u u ) gives 

u M v^(r«") = -v v t (12) 

which can be manipulated to yield 

V (t (ru u )-Vu(Tu (t ) = -Tu (a/ , (13) 

where u^ v = P^P„ (d\u K — d K u\) = is the vorticity of the fluidJU Therefore, for an irrotational 
flow, we can define a potential such that 

Tu^ = d^. (14) 

Thus to solve for irrotational flows of an uncharged fluid, it is sufficient to solve ( I14p and 
ffTTj) . along with an additional equation like fflU]) for every conserved charge. Note that since 
u^u^ = — 1, the equation (1T4"]) may be used to express T in terms of the potential 

r 2 = -(d^Kdy) (is) 

so that (f) determines both and T. 

In general, the charge density q is not related to T . We will, however, restrict ourselves 
to solutions where q/T 3 is a constant. The current conservation equations OH]) are then auto- 
matically solved once the equation (TTTj) is solved. For such flows T(x) is the only independent 
dimensionful quantity that governs the flow. <p{x) determines all the hydrodynamic quantities 
once the ratio q/T 3 is specified. In fact, substituting (fl5|) in ( Tl4|) and finally in (fTTj) one gets 
a single complicated nonlinear differential equation for (j){x^). 

Although the restriction that q ~ T 3 might seem quite ad hoc at this stage, for fluids with 
gravity duals that we will be considering in Section 03), this will turn out to imply that the 
flow is isentropic. Isentropic flows also allow us to parametrically control the temperature of 
the fluid. As discussed above, we need to consider fluids at low temperatures. In the flows we 
consider, derivatives of the velocity, entropy etc. are small, even though their values can and 
should change by 0(1). The equation Qo^|) shows that once we fix the ratio q/s so that T is 
small at some time, it remains parametrically small at all times, since the change of s is of 
order 1. 

9 It can be shown that the condition of vanishing vorticity is identical to the condition d fJl (fu^) = d v {fu^) 
for any scalar function /. In (|14l) / is simply the temperature T . 
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Isentropic sound waves in such a gradient flow are described by small amplitude fluctuations 
of the velocity potential — > 0+50. This induces variations of u M , T and q l , T — > T+5T, — > 
+ 6u^, qi — >■ qi + 5qi. Plugging these into the equations of motion (TH1) . we get 



(T+5T)K + 5<) = c^0 + <9^0 
or u^ST + TSu* 1 = d^5<f). (16) 



Using u^Su^ 1 = we get 



5T = -u^d^l 



Plugging (HZD in (HU) 



and T5w M = P^d v {8(f)) . (17) 

V M [(37^5W + r 3 ^)] =0 
=> ^[v^T 2 (^-2^V)^](50) = O (18) 

For sound waves in a static equilibrium fluid in flat spacetime, this gives (— 3d 2 + d 2 )(6(p) = 0, 
from which we can read off the speed of sound c s = -^=. 

More generally f lTS]) is the Klein-Gordon equation of motion of a massless scalar field in a 
non-trivial background metric, 



yf^GG^dy (50) = (19) 



where 



y/^GG^ = V^gT 2 {g^ - 2u»u u ) 

G^ = V3T 2 [g^ + ^u)j . (20) 

The metric G^ u above, termed the "acoustic metric" , is described by the line element F°l 

ds 2 = V3T 2 |-(1 - ^l 2 )dt 2 - ^p 2 Vidx l dt + (gy + ^j\vj \ dx { dx j ^ , (21) 



where 7(2) = l/yl — v (z) 2 is the Lorentz factor for the moving fluid. If we make the trans- 



formation dr = dt + r 3 a \ dx l , the metric becomes 



ds 



V3T 2 |-(1 - \l 2 )dr 2 + (g^ + -J^ Vi v^J dx'dx^ . (22) 

This is the most general form of the acoustic metric for a relativistic conformal fluid. Note 

that the metric factors vanish or become singular at 7 = which precisely corresponds to 

the speed of sound v = c s = 4^. This indicates that an "acoustic horizon" is formed where 
the flow becomes supersonic and sound waves do not emerge from out of that horizon. 



10 We will use G MJ , and ds 2 for the acoustic metric to distinguish it from the spacetime metric. Here the 
spacetime metric has the form ds 2 = —dt 2 + gijdx l dx^ 
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2.1 Steady flows leading to acoustic horizons 

In this section we will find steady fluid flows with acoustic horizons when the background 
spacetime has a metric of the form 

ds 2 = -dt 2 + dz 2 + R{z) 2 dtt 2 2 . (23) 

R(z) = z corresponds to flat spacetime; we will later consider more general functions R(z) 
and work on spacetimes that are asymptotically flat. We also assume that the thermodynamic 
quantities depend only on z and that v z (z) is the only nonzero component of the velocity. The 
form of the acoustic metric is then 

- \l{z?)dr 2 + d f + R(z) 2 dnl\ (24) 
3 3(l-3 7 (z) 2 ) J 




3 

<£ 7 W 2 (l - V -^)dr 2 + - J;\ t M . + R{zfdn\ 



c 



7 (z) 2 (l - 



The second line can be obtained from the first by using the definition of the Lorentz factor. 
Up to an overall conformal factor the metric is remarkably similar to that of a Schwarzschild 
black hole with a warp factor proportional to (1 — ^(z) 2 ) - an acoustic horizon is present 
at the radius where the flow becomes supersonic. An acoustic Hawking temperature Tjj and 
a surface gravity k can be defined by the standard process of Euclidean continuation of the 
acoustic metric near the horizon; then 



T = — = — 

H 2tt 4tt 



dz 



(25) 



Thermal radiation of quantized phonons is expected from the horizon since Hawking radiation 
is a purely kinematic effect independent of the underlying dynamics |15j . 

In order to get an explicit solution, we need to solve the equations of motion (IT31 and ( TTTT) 



<9,(T 7 ) = T7 = Toc (26) 

d z (R(z) 2 T 3 1 v z )=0 =j> R(z) 2 T 3 1 v z = $ s (27) 

where we have identified the integration constants as the "asymptotic temperature" Too and 
the "entropy flux" $<?. From and fl2ZJ 

^-<) = %J^- (^) 
From the isentropic condition q ~ T 3 it follows that 

q(z) = -f^r • (29) 
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Figure 1: Plot of v z (z) for the spherically symmetric case given by Eqn. (l3TT) with Zh — 1. There 
are two physical branches; neither is valid for z < Zh- The third branch is superluminal. 



2.1.1 Singular radially symmetric solution 

If we take R{z) = z with z G (0, oo), the metric ( )23l) describes flat spacetime. The LHS of 
Eqn. (l28p has a maximum value of for v = c s = 4=, and therefore there is no solution for 
z below a minimum value of 



^=14-5-1 • ( 30 ) 



oo 



3^5- (31) 

This cubic equation for t> 2 has two physical branches, one subsonic and the other supersonic; 
the third branch is superluminal. The physical branches are plotted in Fig. 1. In the subsonic 
branch, as z — > oo, v z — > which is consistent with the identification f l26|) . At the horizon the 
derivatives blow up: 

dv z $ s 2 

-^ = 7f3,3(i-3^)^ 00at z = Zfc - (32) 

As a result, quantities like the Hawking temperature blow up and more importantly the hy- 
drodynamic description breaks down. The pathology can possibly be cured by introducing a 
viscosity [16], i.e. by going to higher orders of the hydrodynamic derivative expansion. 



2.1.2 Solution in more general geometries 

To obtain solutions of first order hydrodynamics that are valid globally, we can choose a more 
general R(z) such that 

• the spacetime is asymptotically flat, R(z) ~ \z\ for z —> ±oo 
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• the maximum value of the RHS of (1281) is ^4=, the same as the maximum possible value 
of the LHS. This condition implies that the minimum value attained by R{z) is 

RmiQ = • (33) 

Then we can construct a smooth solution that changes over from subsonic to supersonic or 
vice-versa every time R{z) attains the above minimum value. For a generic velocity profile, 
the derivatives at the horizon do not blow up because the divergence of at the horizon is 
canceled by the fact that 

dR(z) . . 1 — 3v 1 dv z „ 

-ir = - R{Z) *^W) itz = at z = *■ (34 > 

The simplest example is the wormhole geometry given by 



R(z) = ^ + zl . (35) 

There are two asymptotically flat sheets as z — > ±oo which are connected by a throat at z — 0, 
where R(z) = -R mm = z . From fl33|) we get the condition ^f- = |%, and fl28|) then gives 

^-^ = 72W = iv!^' (36) 

The cubic equation for can be solved for v z . Again there are two physical (subluminal) branches. 
One of them smoothly increases from v z = to v z = 1 for z G (— oo, oo) while the other one 
smoothly decreases. Both solutions have acoustic horizons at z — 0. A plot of the solutions is 
given in Fig. 2. The velocity and its derivatives remain finite near the horizon, as seen from 
the near-horizon expansion of the increasing solution: 

1 V2 z 



V3 ' 3 Zo < 37 > 



and the acoustic Hawking temperature is 

T - 3 



dz 



2\/2tiz q ' 



An obvious problem with this solution is that it reaches the speed of light asymptotically on 
one of the sheets. 

To fix this problem, we can choose, for example 

R(z) = (z 4 -2z 2 z 2 + z A +R 4 mi J (39) 
which is smooth for -R m j n > and has minima at z = ±z where R(z) = R m i n . fl33|) sets 
|^§- = 2 ^| n . This is geometry where two asymptotically flat regions separated by a wormhole 
with two throats. For the velocity profile, we have 

v a _ v *) = 1 = _g R2 ^ uq) 

z{ z) UR{zf 3Vs^-z 2 4 + 4 + RL,' {) 

Now we can find a solution that crosses over from subsonic to supersonic and back to subsonic 
at the horizons and remains subluminal for all z, as shown in Fig. 3. 
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Figure 2: Plot of v z (z) for the wormhole geometry fl35|) with z = 1- There are two physical 
branches, making a subsonic to supersonic (supersonic to subsonic) transition at the acoustic 
horizon at z = 0. 




Figure 3: Plot of v z (z) for the geometry with two throats given by Eqn. (l39l) with z = 1, 
-Rmin — 1- There are two horizons located at Zh = ±^o- The red branch remains subluminal 
for all z. 
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Figure 4: Velocity profile (in red) for the nozzle geometry given by Eqn.( l4"2l) with t> m in = 0.1, 
"Wmax = 0.9, Zq — 1 and the corresponding R(z)/R min (in green). The minima in R(z) at 
z = ±z correspond to v — c s . 



2.1.3 Nozzle geometry 

We can also choose the spatial section of the metric to be asymptotically cylindrical with 
R(z) — > constant at z — > ±oo. For variety we consider geometries with a toroidal cross-section 

ds 2 = -dt 2 + dz 2 + R(z) 2 (d6l + de 2 2 ). ( 41 ) 

The form of the acoustic metric ( 123]) remains very similar with the replacement of dQ 2 , by 
d8f + dQ\. The acoustic Hawking temperature is still given by ( J25l) . As an example, we can 
take the following profile for v (z) and then solve (1281) for R(z): 

v z = v min + (w max - v min ) sech ( — J . (42) 

The f m ; n term ensures that R{z) — > constant and remains finite for large \z\. The profiles of 
v (z) and R{z) are shown in Fig. 4. There are two horizons located at 

z± = ±Zq cosh . [AS) 

C-s ^min 

Assuming v z > 0, the fluid passes the speed of sound at the left horizon Z- and again returns 
to subsonic speeds as it crosses the right horizon z + . 

2.2 Sound waves 

In this subsection we examine the behavior of sound waves around the background flows 
described in the previous subsections. Sound waves are fluctuations of the velocity potential 
cf) which satisfy a massless Klein-Gordon equation in the background acoustic metric (124"1) . We 
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will consider background flows on a simple wormhole geometry - more complicated wormholes 
or nozzle geometries can be treated along similar lines. 

Near the acoustic horizon (chosen at z — 0), the metric is similar to the usual Schwarzschild 
metric, so we expect a large blueshift effect for outgoing modes. To display this, it is sufficient 
to consider modes in the s-wave. Let us first rewrite the metric in terms of null coordinates 
u, v as follows 

d~ s 2 = V3T^(z) - ^ 2 )dudv + R\z)(d6 2 + sin 2 #d0 2 ) J , (44) 



where 



dz \rfvdz dz dz 
du = dr = = dt + - — s— ; = = at 



v^l-fT 2 ) l-§7 2 v/3(l - |T 2 ) *>+(*) 



dz , -^fvdz dz , dz , iK . 
dv = dr+ —= = dt + k-t + = = dt — . (45) 

Here v± (z) = i±v(z)c * s ^ ne relativistic sum of the local fluid velocity and the velocity of sound 
c s — 1/ v3- Then close to the acoustic horizon R(z) — > R(Q) is finite, and the s-wave solutions 
of the Klein-Gordon equation are approximately 

„ e ~ iulu ifj_ ~ e~ iujv . (46) 

In the asymptotic region, where the velocity becomes constant, we obtain the usual spherical 
Bessel functions. To analyze the near-horizon behavior, we can expand the velocity field near 
the horizon as v(z) ~ c s + |/cz + . . ., where k = 1 1 is the surface gravity at the horizon. 
Using this in (146]) . we find 

^ + ~ e -^-^ + ° {z2) ] (47) 
^_ ~ e -H t -^ + 75 + °^} (48) 

The ip + mode is continuous at the horizon and is right moving with a velocity which is the 
relativistic sum of ^= with itself. The ip- mode has rapid oscillations near the horizon: 

z<0 : ^_ ~ e -M*-£M-*)] (left-moving) 
z>0 : ^_ ~ e -^[*-^ ln ( 2 )] (right-moving) 

indicating that inside the horizon, both modes are right-moving. 

To extend these modes away from the horizon, we employ an eikonal approximation. We 
decompose the sonic fluctuation into a rapidly varying phase or "eikonal" (here A ^> 1) times 
a slowly varying envelope 

^) = A(x Ai )e- iA5( ^ ) (49) 
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and plug it into the wave equation <9 M [\/-GG F 9 1/ ] ip = with the acoustic metric (j2"Tl) written 
out in the original t-z coordinates 

d§ 2 = r 2 (-(1 - ^7 2 )rft 2 - ^-fv z dtdz + (1 + f-y^cfe 2 + R{z) 2 {dB 2 + sin 2 ^ 2 )) . (50) 

We then get a sequence of differential equations for Six' 1 ) and y4(x M ) by expanding the wave 
equation order by order in A: 

0(A 2 ) : d^S{x a )d^S{x a ) = (51) 
0(A) : 2d ll S{x a )d^A{x a ) + A{x a )V 2 S{x a ) = (52) 



The leading equation (15T|) can be used to solve for S(x a ). Let us consider s-wave solutions 
independent of 9, 4>. With the ansatz 

XS{x a ) =ut- f(z) (53) 

we can solve for the phase, 

AS±(x a ) = wt + w / dz- T -^- = ut-u dz ^ (54) 

J I-37 2 7 

The momenta of the wavepackets are given by derivatives of the eikonal["1 p,, = — d^[\S(x a )} 
giving 

p t = u 

1±vc s w , , 

P 2 = w — — — = — , (55) 
V ±c s v± 

where v± has been defined above. The upper and the lower signs would normally correspond 
to right- and left-moving sound modes. However, if the fluid (assumed to be right moving) has 
a velocity greater than the speed of sound, then both the modes become right-moving. Using 
the ansatz A(x a ) = A(z) in the subleading equation (!5"2i) . a full solution is obtained: 

ip± = e L J (56) 



3 Gravity dual of acoustic solution 

The fluid-gravity correspondence [31 EH El El CO EH [22] provides a correspondence between 
solutions of certain fluids with classical solutions of suitable Einstein-Maxwell equations in a 

4 + 1 dimensional spacetime with a negative cosmo logical constant. In our case the fluid is a 

n The leading order equation (|51l) is thus a null geodesic equation, = for the phonons. 
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conformal fluid with a global U(l) charge and the higher dimensional bulk theory is given by 
the five-dimensional action 



S = / d b x J=B 



R + 12- F AB F AB - ^e EABCD A E F AB F CD 



(57) 



where G is the five dimensional Newton constant, the indices A,B run from to 4, A B is a 
U(l) gauge field, and we have chosen units in which the cosmological constant is A = —6. The 
above action is a consistent truncation of IIB supergravity for k = 1/(2a/3). We will, however, 
allow arbitrary values of n. 

A uniform charged black brane solution of this action is given in a boosted reference frame 

by 

ds 2 = —2u^dx^dr — r 2 V(r,m,q)u fJj u u dx IJ 'dx u + r 2 P lll ,dx tJ 'dx L ' 
A = T^u,dx» (58) 

where w M are constant 4-velocities (the indices /i, v = • ■ ■ 3) and P^ = r/^ + u^Uy is the 
spatial projection operator. The function V(r, m,q) is given by 

~2 

771 Q 

V(r,m,q) = l-- + q - (59) 

where m and q are parameters of the solution and we are using the notation of [22J. 

This solution is dual to a charged fluid in equilibrium living on the flat boundary of the 
five- dimensional spacetime. The fluid is strongly-coupled M = 4 SU(A^) Yang-Mills theory, 
viewed in a boosted frame with coordinates The temperature T, charge density q, energy 
density e and entropy density s of the fluid are given by [22J 

R ( q 2 \ 1 
T = — - ( 2 — — r , q = V3aq, e = 3am, s = AnaR%, a = — — (60) 

27T \ J 167TG 

where R + denotes the radius of the outer horizon, i.e. the largest root of the equation 
V(r, TTi, q) = 0. The energy momentum tensor T^ u and the charge current J M of the fluid 
are given by 

T ^ = + 4m a^) J v = V u v ( 61 ) 

The expressions (!60|) and (IBTj) involve the bulk parameter G. In our units, G is related to the 
rank of the gauge group of the boundary theory by 

G= 2^> a= ltG = ^ (62) 
With the substitutions in ( pDj) . the equation of state e(s, q) becomes identical to the condition 
V (R + (s) , m(e) , q(q)) = 0. The equation of state for a charged conformal fluid with a gravity 
dual is thus 

^^r + sJw^ ^ 1 *^^^)* (63 > 



s \ 3 q ( Ana 
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The temperature and chemical potential can be obtained by taking derivatives of e(s, q) using 



7r \4TtaJ \ 3s 2 

This temperature reproduces the value quoted in fl60l) . In the uncharged limit i? + = 7rT, and 
one can fix the value of c defined in (J7J) by comparing it with (I6"3"|) and requiring that T = T 
for uncharged fluids; this gives c = a7r 4 and hence 

6^=0 = 3a(7rT) 4 , s M= o = 47ra(7rT) 3 . (66) 

For charged fluids at finite \i there is a zero-temperature limit, reached when R + = ^/(2y6) 
and T = /i/(192*7r). In this limit, 

<Zt=o = ^ , e T=0 = ^ , s T = = ^ . (67) 

As discussed in Section 2 we restrict our attention to isentropic flows. For such flows q/T 3 
is constant and since e = 3«7r 4 T 4 , it follows from the equation of state (163!) that T/R+ is a 
constant, and therefore that for such flows the entropy per unit charge s/q is constant. The 
first equation in (1651) shows that by choosing 

1 -^(£) 2 « 1 

we can keep the temperature T -C -R+ everywhere and at all times. 

The gravity dual of a general fluid motion is then constructed in a derivative expansion 
as follows. First, we replace the parameters of the solution by functions of the boundary 
coordinates x M , w M — > u^(x),m — > m(x),q — > q(x) which respectively represent the velocity 
field, energy density field and the charge density field of the fluid. We also replace the flat 
boundary metric rj^ with a curved metric g^ u (x). With these replacements, ( )58|) is no longer 
a solution of the bulk equations of motion. Second we need to add correction terms to the 
metric and the gauge field so that the full metric and the gauge field now solve the equations 
of motion. This second step is of course impossible to perform in an exact fashion. However, 
these corrections can be calculated systematically in a derivative expansion, provided that the 
derivatives of m^(x), m(x), q(x) with respect to x v are small compared to the outer horizon 
radius R + . To lowest nontrivial order in the derivative expansion, the modified metric and 
gauge fields are 

ds 2 = —2u^dx^dr — r 2 V \r,m,q)u il u v dx >l dx 1 ' ' + r 2 P ilv dx p 'dx v 

2 2r 2 

+ -r(V a u Q )u fl u v dx' 1 dx v + —a fMU F 2 (p, M)dx»dx v 

3 R^ 

- 2ru IJ u a {y a u u )dx^dx v 

2u, ( + ^(P u x d x q + 3(u A V A u,,)g)*i(p, M) ) dx»dx v , (69) 

\ mr 4 R' + J 
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A 



where we have defined the quantities 



dx» (70) 



M =i Q =n p= k (71) 

V M is a covariant derivative with boundary metric g^, and 

^ = 9^e va ^u u V a up , a, v = Ip^p-P (V a u p + V p u a ) - ^ u (V a u a ) . (72) 
The functions Fi(p, M, Q) and F 2 (p, M) are defined as 

Fl(p , M , Q) = 1 (i _ ^) ^ p _^ ( 1 _ ^ (1 + J_)) , (73) 

Note that in the above expressions m, q, R + , M,Q, p are also functions of the boundary coor- 
dinates x^ since m and q are functions of x M . 

This is a solution of the bulk equations of motion, provided that m(x), q(x) and u^x) are 
such that the energy momentum tensor and current 

_, mix) , , . . . N , T VSqix) , x /__x 

^ = IBttG^^ + 4w a*W w,/ W) J <" = i67rG ^ 

are covariantly conserved, 

V^T"" = V M J M = (76) 
Thus every solution of fluid dynamics leads to a bulk solution. 



3.1 Gravity duals of dumb holes 

We now apply the results of the preceding subsection to construct gravitational duals of the 
fluid flows with acoustic horizons that were studied in Section 2. These flows are special in 
several ways: first, the background spacetime metric of the fluid is of the form (T23]) or ( I4ip 
where the only inhomogeneity is in the z direction. Second, both the background flow and 
the sound wave fluctuations have vanishing vorticity. Third, the background flows as well 
perturbations around them are isentropic. 

It follows from the isentropic condition that the quantities M and Q are constants. As ar- 
gued above (see discussion following equation ({15]) ). for isentropic flows there is just one length 
scale, and all quantities are related to this length scale by dimensional analysis. In particular, 
the dimensionless quantities M and Q must be constant. In addition, the inhomogeneous parts 
of all quantities which appear in the bulk metric and the gauge field are determined in terms 
of a single scalar field 4>(x). 
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As discussed in the introduction, in order for the acoustic Hawking radiation to be de- 
tectable we need to consider fluids which have a very small ambient temperature. This means 
that the constant quantities Q and M need to be close to their extremal values, Q ~ \/2 and 
therefore M ~ 3. While the various quantities like e(x),q(x) can change by 0(1) amounts 
(only their derivatives are small), the isentropic condition ensures that if the fluid temperature 
is initially small it will remain small (see the discussion in Section 2 above). 

To construct the background fluid flow, we simply need to insert the velocity potential 
0o for the solutions of Section 2 into the general metric and gauge field in (159"]) and (170]) . 
The conditions of vanishing vorticity and isentropic flow simplify these general expressions 
somewhat. The most drastic simplification appears in the expression for the gauge field, 
equation flTOj) . In fact the first order corrections (in the derivative expansion) to the gauge 
field vanish for isentropic gradient flows. To see this we note first that V ' a up can be replaced 
by daUp in the expression 1^ of f l72|) . Then using Tu^ = d^cf) we get 

J„ = g^ a ^[-^d a Td^ + ^d a d^] = (77) 

due to antisymmetry of the epsilon symbol. The third term on the RHS of (170]) also vanishes, 
as can be seen by applying the isentropic condition q/T 3 =constant and the relations (|T4|) and 
(1151) to the expression 



Pft x q + 3(u x \7 x u^q = 3aT [Td^T + Tu x V X {T 

1 



3aT 



2 



d^d a <pd a <p) + d x <pv x d, 



(78) 



Thus, to first order in the derivative expansion, the bulk gauge field is given by the first term 
of the right hand side of (1701) . which is just the term which would have resulted from a simple 
boost of the original black brane solution. In our case the charge density and the 4-velocity 
appearing in (170"]) are functions of z, as determined by the fluid flow on the boundary. So there 
is a nonzero electric field component along the z direction, given by 

F 0z = - ^°° v z d z v z (79) 

where we have used ( 129]) to express q(z) in terms of the velocity v z . 

The expression for the bulk metric simplifies as well. In the fourth line of ( 169]) . the first 
term is proportional to Z M which vanishes for our flows. The second term is proportional to 
defined in (175]) and vanishes as shown above. 

In the derivative expansion, the relationship between the boundary and the bulk becomes 
essentially local. The bulk solution can in fact be constructed approximately by patching 
together tube geometries obtained by extending the boundary data in a given region of the 
boundary to the bulk using the radial equations of motion. Consequently we expect that the 
acoustic horizon of the fluid flow on the boundary extends trivially into the bulk. We will 
explicitly verify this in the next subsection. 

However, this tubewise approximation breaks down in regions where the local geometry is 
exactly extremal. This is apparent in the results of [21] and (22]. Furthermore, recent work 
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on perturbations around extremal black holes shows that the relevant low energy expansion is 
different from the naive derivative expansion [23]. Our results hold close to extremality, but 
not exactly at extremality. 



3.2 Gravity duals of phonons 

The gravity duals of phonons in the fluid are quasinormal modes of metric and gauge field 
perturbations. Once again, construction of these modes is trivial. We need to write 



compute T(x) and hence m(x),q(x),u fl (x) in terms of 5<p, substitute into and (170]) . and 
consider the terms which are linear in (3. By construction, these modes satisfy ingoing boundary 
conditions at the bulk horizon. 

The fluctuations of the gauge field obtained by this procedure have a particularly simple 
form 



For the background flows considered in Section 2, we have found solutions to the wave equation 
(fl"8l for 5(j) in the region close to the acoustic horizon. Upon inserting these solutions into (IHTj) . 
we see that the fluctuations SA^ have a characteristic behavior near the acoustic horizon, viz. 
ingoing waves are smooth while outgoing waves have rapid fluctuations. This is the precise 
sense in which the fluctuations perceive the acoustic horizon, which has now extended into the 
bulk. From the nature of the solution that the extension of the acoustic horizon into the bulk 
is rather trivial - i.e. the horizon perceived by these modes is at the same value of z as the 
acoustic horizon on the boundary, and for all values of r. 

The fluctuations for the components of the metric can be similarly worked out and also 
see a horizon structure at the same value of z. We therefore conclude that there are certain 
quasinormal modes of the bulk metric and the gauge field which perceive a horizon. If these 
bulk modes are quantized, one should find a thermal bath of such modes characterized by the 
temperature of the acoustic horizon on the boundary. 

4 Regime of validity 

It is important to check that the fluid flow described above is consistent with the standard 
conditions for validity of hydrodynamics. Roughly speaking, hydrodynamics is valid when the 
gradients of velocities, temperature and charge densities are small compared to the inverse 
mean free path l m . For charged conformal fluids considered above, there are two scales - the 
temperature T and the chemical potential /i = z/T, so that l m ~ f{y)jT. The function f{y) 
is of order one for generic values of z/, but there is an upper bound on z/, v c where f{y) has a 
simple zero. It is possible to take the limit of v — > v c simultaneously with T — > such that 
l m is finite - the dual of this is in fact the extremal black hole. For the flow described in the 




(80) 



5A, = -pl—lidMi&^S* + (d^oX^o)] W) 



(81) 
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previous section, in this limit we have 



dv z (z) 



dz 
1 dT{z) 



T dz 



1 

<<: r 

'rn 

« r (82) 



In particular, since the acoustic Hawking temperature Tjj is ^ \^~\ z > this implies that 

Th y- (83) 

For observable acoustic Hawking radiation, the Hawking temperature should be higher than 
the fluid temperature. So we require 

T<T H (84) 

Furthermore, the frequency of the sound waves should also be small compared to the basic 
scale, u <C l// m . However, the finite Hawking temperature is going to introduce an upper 
bound on the allowed wavelengths, due to periodicity in Euclidean time; thus u > T H . Thus 
we need 

T < T H < to < (85) 

For fluids with no conserved charge, there is only one energy scale, namely, the temperature 
T; thus l/l m ~ T and ( I84p cannot be satisfied. Although the solution is otherwise valid, the 
Hawking radiation, at a temperature much lower than the ambient temperature, is not going 
to be observable. For fluids with a conserved charge the condition (1841 does not pose a problem 
because now we have two length scales, the temperature T and the chemical potential /i. For 
fluids very close to zero temperature, the mean free path will be governed only by fi. We can 
thus have 

0^T<T H <u^^-^T^fi (86) 

I'm 

The ability to construct a gravity dual using a derivative expansion imposes further con- 
ditions. In the presence of a nonconstant R(z), the validity of the derivative expansion of 
the solutions of the bulk equations of motion requires that the curvature invariants and all 
invariants constructed out of the derivatives of curvature be small. An example of such an 
invariant is g^V \JZV ' U TZ, where TZ is the Ricci scalar and we require for this example 

(<rV^V„7^ < p. (87) 

I'm 

We get additional conditions if some of the boundary directions are compactified as in the 
nozzle geometry of Section (I2.1.3p . If one boundary direction of a AdS x S geometry is made 
compact with a radius R(z), the dual is an AdS soliton |25j which caps off the geometry at 
a value of the radial coordinate r = l/(2R). For a black brane geometry, compactification of 
a boundary direction would lead to a similar modification of the usual black brane geometry. 
However if R(z) ^> l/(2i? + ), where R + is the location of the black brane horizon, the place 
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where the bulk geometry would cap off is far inside the black brane horizon. In this situation we 
can continue to use the standard black brane geometry with a compact longitudinal direction. 
We will therefore require that for all z, 

R{z) > 1/R+ (88) 

for the solution of Section (I2.1.3p . Finally we require for the nozzle solution that R(z) be finite 
for large z. The geometry is then asymptotically R x T 2 and has an AdS dual. 



4.1 Validity of our solutions 

Finally, we determine the range of parameters for which our approximations are valid, for the 
specific flows studied in Section 2. Let us first discuss the wormhole solution of equations ( 1391) 
and PU|) . The solution has four parameters, %o, goo, zq and -R mm . $s is fixed by (1331 once -R m i n 
is chosen. Since v < 1, 7 remains finite for all z. T = — > and we can have a valid derivative 
expansion w.r.t. To- In order that the curvatures are small, we require R m i Q 3> l/To- The 
derivatives ^7, 7=^7 and the derivatives of the curvature are all proportional to i. Thus we 
require ^ <C %o- Qoo can be chosen such that we are always at very low temperatures, following 
the discussion around equation (|68|) . To summarize, the conditions for our wormhole solution 
to be valid are: 

T^O, -<w<7^, -^<Too. (89) 

Zq iX m i n 

For the nozzle solution described by (jl2l . the parameters are z , f max , f m m, Tx>, qoo and $5. 
As in the previous case, all derivatives in the solution are proportional to We need t> max < 1 
for the derivative expansion to be valid (so that 7 remains finite and T remains non-zero) and 
we obtain the same conditions as in fl89|) . The condition fl88|) is same as the requirement that 
the background curvature remains small: -R m i n ^> l/Tx,- Since -R m i n is given in terms of Too 
and $5 by fl33|) . this implies $5 ^> Too- Moreover, we need v m i Q > in order that R(z) be 
finite at large z - the asymptotic geometry remains Kx T 2 , and we have an asymptotically 
AdS gravity dual. Summarizing, the conditions for validity of our nozzle solution are: 

T ->• 0, — < u < To < $5, < v min < w max < 1. (90) 
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